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ABSTRACT 



1^' Context. Young stellar objects (YSOs) occasionally experience enhanced accretion events, the nature of which is still poorly under- 

Q I stood. The discovery of various embedded outbursting YSOs has recently questioned the classical definition of EXors and FUors. 

^ , Aims. We present a detailed spectroscopic investigation of the young eruptive star PV Cep, to improve our understanding of its nature 

j/3 ■ and characterise its circumstellar environment after its last outburst in 2004. 

C^ ' Methods. The analysis of our medium-resolution spectroscopy in the near-lR (0.9-2.35 yum), collected in 2012 at the Large Binocular 

Telescope with the IR spectrograph LUCIFER, allows us to infer the main stellar parameters (visual extinction, accretion luminosity, 
mass accretion and ejection rates), and model the inner disc, jet, and wind. 

Results. The NIR spectrum displays several strong emission lines associated with accretion/ejection activity and circumstellar envi- 
1^ ' ronment. Our analysis shows that the brightness of PV Cep is fading, as well as the mass accretion rate (2x 10"^ Mq yr"' in 2012 vs 

l ' ~5x 10"^ Mq yr"' in 2004), which is more than one order of magnitude lower than in the outburst phase. 

Among the several emission lines, only the [Fe ii] intensity increased after the outburst. The observed [Fe n] emission delineates 
blue- and red-shifted lobes, both with high- and low-velocity components, which trace an asymmetric jet and wind, respectively. 
The observed emission in the jet has a dynamical age of 7-8 years, indicating that it was produced during the last outburst. The 
^h I visual extinction decreases moving from the red-shifted (Av(''e(i)=10. l±0.7mag) to the blue-shifted lobe (Av(WMe)=6.5±0.4mag). 

f~^ , We measure an average electron temperature of 17 500 K and electron densities of 30 000 cm""* and 15 000 cm"'' for the blue and the 

rf^ ■ red lobe, respectively. The mass ejection rate in both lobes is ~ 1.5x10"^ Mq yr"', approximately matching the high accretion rate 

observed during and immediately after the outburst (Mo,„/M„^.j.~0.05-0. 1). The observed jet/outflow asymmetries are consistent with 
an inhomogeneous medium. 

Our modelling of the CO emission hints at a small-scale gaseous disc ring, extending from ~0.2-0.4 AU to ~3 AU from the source, 
isj^ . with an inner temperature of ~ 3000 K. Our H i lines modelling indicates that most of the observed emission comes from an expanding 

rS ' disc wind at re=10000 K. The line profiles are strongly affected by scattering, disc screening, and outflow self-absorption. 

Conclusions. According to the classical definition, PV Cep is not an EXor object, because it is more massive and younger than typical 
EXors. Nevertheless, its spectrum shows the signature of an 'EXor-like' outburst, suggesting a common origin. 

Keywords, stars: formation -stars:circumstellar matter -stars: pre -main sequence -lSM:individual objects: PVCep - Infrared: ISM 

1. Introduction 



Send offprint requests to: A. Caratti o Garatti, e-mail: 

acaratti@mpifr-bonn.mpg . de Young stellar objects (YSOs) are characterised by accretion and 

* Based on observations collected at LBT. The LBT is an interna- ejection processes, which are closely related to and accompany 

tional collaboration among institutions in the United States, Italy and the whole star-formation period from the protostellar to the pre- 

Germany. LBT Corporation partners are: LBT Beteiligungsgesellschaft, ^ain sequence phase, with highly variable strength and duration. 
Germany, representing the Max-Planck Society, the Astrophysical 

Institute Potsdam, and Heidelberg University; The University of 

Arizona on behalf of the Arizona university system; Istituto Nazionale Coiporation, on behalf of The University of Notre Dame, University 

di Astrofisica, Italy; The Ohio State University, and The Research of Minnesota and University of Virginia. 
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In addition to continuous accretion, YSOs may occasionally 
experience enhanced accretion events, which produce episodic 
increments of their optical and infrared brightness. These ob- 
jects are called young eruptive stars, usually divided into FUors 
(named after FU Orionis, the prototype) and EXors (from EX 
Lupi). FUors exhibit a large brightening u p to ~5-6 mag in the 
optical lasting from few years to decades (Hartmann & KenvonI 
L1996), whereas EXors (Herbig 1989) display smaller bursts (1- 
3 mag) with shorter duration (from days to months), and higher 
frequency (every few years). The two subclasses also differ 
from each other in the observed spectral features in the opti- 
cal/NIR regime. Namely, EXors have emission-line spectra, sim- 
ilar to those from Classical TTauri stars (see e. g., Herbig 2007), 
whereas FUors show absorption-line spectra, which render dis- 
tinct spectral types (SpTs) as a function of the wavelength (i. e. 
F or G SpT s can be inferred in the o ptical, M SpT in the 
NIR; see e. g.. lHartmann & Kenv on 1996). Moreover, EXors and 
FUors have been historically classified as low-mass, pre-main 
sequence objects. However, the discovery of various embedded 
outbursting YSOs, especially in the NIR (e.g., OO Ser, V1647 
Ori, HBC 722, V 2492 Cyg , V2775 Ori; see e. g. [H odapp et al. 
1996; Fedele et al. 2007; Kospal et al. 201 lat lCm^ et alJ201 U 
Caratti o Garatti et al...201 1,) . suggests that these processes may 



take place in YSOs within a wide range of masses and ages and 
that there might not be a neat separation between the two sub- 
groups. 

These outburst events are likely produced by disc instabili- 
ties, and YSOs may accrete a sign ificant ammount of their mass 
during such episodes (see, e. g., lEvans et al] 120091; IVorobvovl 
1200 9*). In addition, enhanced accretion increases the outflow 



activity (Cabritetal. 1990; Britt ain et al. ' 2007). The material 
ejected by such outbursts might produce bright knots along 
the jet by tr iggering MHD instabilities in the flow (see e. g., 
iFendll l2009l) . Like the majority of YSOs, also many young 
eruptive stars drive outflows, HH or jets. For example, V2492 
Cyg and V346 Nor are known to drive Herbig-Haro ob- 
jects (HH 569 and HH 57, respec tively; see iBallv & ReipurM 
l2003l; lReipurth & Wamstekerfll983l) . On the other hand, OO Ser 
is ass ocia ted with o utflow activity and shock-excited H2 emis- 
sion jHodapp et al] |2012|) . whereas Parsamia n 21 shows Ha 
knots with dynamical ag es of 40-80 years ( Staude & Neckell 
Il992t lK6spaletaI]l2008l) . However, it is diflicult to associate 
specific eruptive/outburst events with knots or augmented ejec- 
tion along the flows. To the best of our knowledge, Z CMa is the 
only case where ejected knots along the flow could be associated 
with observed outbursts (IWhelan et al.ll2010 ). Thus, the study of 
such phenomena is of extreme interest, providing us with infor- 
mation about the accretion/ejection mechanism as a whole. 

Among the outbursting sources, PV Cephei is a unique ob- 
ject. It is a pre-main sequence young eruptive star, strongly 
variable, located in between the L 1 155 and L 1 158 clouds . 
Its distance is debated, betwee n 325 pc (IStraizvs et al.l Il992h 
and 500 pc (iCohen et al.l ll98lL the latter value is assumed 
throughout the paper), as weU as its spectral type. Although 
largely recognised as an A5 (because of the Hi Unes observed 
in absorption; see e .g., ICohen et all ll98U iThe et all [19941 
[Abraham et al.l I2OOU) . an absorption spectrum of Sp T G8-K0 
has also been reported (iMagakian & Movsesianl200lh. The esti- 
mated mass of the central source ranges from 2.4 Mq dKun et alJ 
l2011h to ~4Mq (Hamidouche 2010) with a massive accre- 
tion disc of ~O.8M0 (Hamidouche 2010), which drives a pre- 
cessing jet with Herbig-Haro objects (HH215, HH315 and 
HH415; Reipurth et al. 1997) and a massive asymmetric CO 
outflow (Arce & Goodmanii2002bi) . Adaptive optics (AO) ob- 



servations exclude the presence of a close companion (down 
to -50 AU; IConneUev et al.ll2009 1. which could trigger the out- 
bursts. 

PV Ce p was spectro scopically and photometrically classified 
as EXor (lHerbig||1989l) . after the outburst in 1977. After ex- 
hibiting a large outburst in 2004, its brightness faded by sev- 
eral magnitude s (up to four in the I and R bands) between 
2005 and 2009 jKun et alJ |2( )Tlh . with a transient peak in 2008. 
Its photometric decline was caused by the r educed accretion 
rate and the increased circumste llar extinction (Kun et al.ll201 ll ; 
iLorenzetti et al.l201ll) . Recently. lLorenzetti et al.l(l201 ll) pointed 
out that PV Cep is not a genuine EXor, being more massive and 
complex than typical solar-mass EXors. 

To clarify the nature of PVCep, we present its medium- 
resolution spectroscopy in the near-IR (0.9-2.35 //m), collected 
at the Large Binocular Telescope with the infrared spectrograph 
LUCIFER. 

This paper is organised as follows. Section |2] illustrates our 
observations and data reduction. In Sect. |3] we report the results 
with a detailed description of the detected spectral features, and, 
in Sect.|4] we derive the main accretion and ejection properties 
of the YSO and disc physical parameters from the analysis of its 
[Fen], Hi, and CO features. In Section |5] we then discuss the 
accretion/ejection properties in PV Cep, its peculiar asymmetric 
jet and its circumstellar environment. Finally, our conclusions 
are drawn in Sect. |6] 



2. Observations and data reduction 

Our spectroscopic observations of PVCep were acquired at 
the Large Binocular Telescope (LBT) on Mount Graham 
International Observatory (AZ, USA), using the infr ared camera 
and spectrograph LUCIFER at medium resolution (iSeifert et al.l 
.2003.) . The N1.8 camera was used, corresponding to a spatial 
scale of O'.'25/pixel. The same spectral dataset was taken on 
the 23rd and 24th of June 2012 (with an average seeing of 1" 
and 0'.'8 seeing, respectively) and covers the z (0.88-1.01 //m), J 
(1.16-1.31jum), H (1.54-1.74yum) and K (2.01-2.35yum) spec- 
tral segments, with a total integration time (per night and per seg- 
ment) of 1080, 1680, 480 and 192 s in the z, J, H and K bands. 
We adopted the 210_zJHK grating unit with a 0'.'5 slit width, 
leading to a nominal spectral resolution of K ai6900, 8500, 
7800 and 6700 in the z-, J-, H- and /T-bands, respectively. The 
slit was positioned on the target along the main outer outflow 
axis (P.A.=348°). All the raw data were reduced by using the 
IRAlQ package and applying standard procedures. Each obser- 
vation was flat fielded, sky subtracted, and corrected for the cur- 
vature derived by long-slit spectroscopy, while atmospheric fea- 
tures were removed by dividing the spectra by a telluric standard 
star (HD 207636, AO spectral type), normalised to the black- 
body function at the stellar temperature, and corrected for its 
intrinsic photospheric absorption features. The raw spectra were 
wavele ngth calibrated using t he bright OH lines detected on each 
frame JRousselot et aLlBOOOl) . with average uncertainties of 0.2- 
0.4 A, depending on the number and intensities of the OH lines 
detected in the considered spectral segment. The average val- 
ues of the instrumental profile in the dispersion direction, mea- 
sured from Gaussian fits to the OH sky lines, is 2.1, 2.2, 2.9, 
4.9 A, in the z, J, H, and K band segments, respectively. Radial 



' IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated 
by AURA, Inc., cooperative agreement with the National Science 
Foundation. 
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Wavelength (A) 



ible on the continuum. For each detected line, Table [T| (shown 
in Appendix) reports the observed (vacuum) wavelength (Aohs), 
the calibrated flux (uncorrected for the extinction) and its un- 
certainty {F,AF), the equivalent width (W,0, the full width half 
maximum {FWHM), the observed signal to noise ratio (S/N), 
the line identification (ID), its vacuum wavelength (A^ac), and 
the resulting radial velocity (Oi (peak)). The equivalent widths 
and line fluxes were calculated by integrating across the line, af- 
ter subtracting the continuum. A few lines could not be properly 
identified and are labelled as uncertain (ID + '?') or unidentified 



24000 ('?') in Tab.rri 



Fig. 1. LUCIFER z-, J-, H- and K-band medium resolution flux- 
calibrated spectrum of PV Cep. 



velocities of the observed lines were also obtained through sin- 
gle or multiple Gaussian fits in case of blended lines, measured 
in the observer local standard of rest (LSR) and corrected for 
the cloud speed with respect to the LSR (udoud = -3.0 km s"', 
[Torrelles et al. 1986). As a result, the velocity calibration is ac- 
curate up to ~10-20kms ', depending on the considered spec- 
tral segment. Finally, each spectral segment was flux calibrated 
through a spectrophotometric standard star (Hip 106737) ob- 
served using the same instrumental settings and slit widths of 
O'.'S and I'.'S, to estimate and correct for the flux losses due to the 
seeing. 



The most prominent lines in the PV Cep spectrum are cir- 
cumstellar features, originating from YSO ac cretion or inner 
winds, such as, e. g.. Hi, Can (see e.g ., Muzer oUe et al.|[l998bl: 
iNatta et al.l l2004t ICaratti o Garatti et al.. ,20121) . s hocks along 
the jet, such as, e.g., TFeiil an d H2 (see e.g., iNisini et al] 
I2OO2I: ICaratti o Garatti et al."2006'), inner disc region (e. g.,N ai, 
CO) (see e. g.. Kospal et al. 201 lb: Lorenzetti et al] [201 ih . or 
chromospheric activity (e .g., Fei, Fe n, Mg i, Ci) (see e.g., 
iHamann & Perssonlll992allb[lKellv et alJll994l) . 



A few fluorescent emission lines from Fen (i.e., at 1.0 
and 1.69/im) and Ci (i.e., 1.17-1.19/im) have also been 
detected. They are usually pumped by UV photons, likely 
from the Lya or continuum emissio n of the stellar photo- 
sphere (see e. g., iBautista et al.' '2004^ 'Johan sson & Letokhovl 
l2004tlLumsden et al.ll2012n Wal mslev et al. 200^^ 



3. Results 

First, the two spectroscopic datasets were analysed separately to 
check for variability. There is no clear evidence of variability in 
both line profiles and YSO continuum between the two nights, 
within the uncertainties (AF/F) given by the flux calibration of 
our spectrophotometric standard star (about 20% in the / band, 
and between 5% and 10% in the J, H, and K bands) Therefore, 
we combined the two datasets to increase the signal to noise ra- 
tio of our spectrum. The final NIR flux-calibrated spectrum of 
PV Cep is presented in Figure [1] and shows a steeply rising con- 
tinuum with several strong emission lines. The flux-calibrated 
spectrum translates into the following /, J, H, and K band mag- 
nitudes: 14.7±0.3, 12.30+0.08, 10.73±0.08mag, and8.80±0.04, 
respectively. These values are smaller than those previously ob- 
served before, during, and afte r the outburst in 2004 or the tran- 
sient peak in 2008 (see e. g., iLorenzetti et al.l 1201 ll: iKun et all 
I2OIII) . indicating that PV Cep's brightness is slowly fading, es- 
pecially in the J, H, and K bands (Am/~1.5 mag, Amy~2.1 mag, 
Am/f ~2.6 mag, and Am^~2.3 mag, since 2008). 

3. 1 . Detected lines 

Figure |2] displays our PVCep continuum-normalised spectrum, 
labelling the most prominent features. Lines were identified 
using wavelengths from the NIST Atomic Spectra Database 
Linefl In addition, NIR spectra from se veral active Classical 
TTauri stars (CTTSs) and EXors (e.g.. iKospalet al.1 l201Ibl: 



iGarcia Lopez et al.l20I itlCaratti o Garatti et al.l20I2l) were also 
scrutinised to recognise the observed features. We detect sev- 
eral forbidden and permitted lines from atomic and molecular 
species, but no obvious absorption photospheric feature is vis- 



The detected emission lines are normally observed in the 
NIR spectra of EXors and very active CTTSs, although the fluo- 
rescent lines are more typical of Herbig Ae/Be stars. Some of the 
brightest lin es, namely H i, Mg i, Na i, CO bandheads, were also 
detected by ILorenzetti et al.l (l2009l) during the transient peak 
in 2008 through NIR low-resolution CR a!250) spectroscopy. 
As for the continuum, absolute line fluxes also show a signifi- 
cant decrease of about one order of magnitude, with the excep- 
tion of the [Fe n] emission lines, which were barely vis ible in 
the NIR low-resolution spectrum (K a: 1200) of Connelle v et al.l 
(2009), taken in July 2007, a nd below the d etection thresh- 
old (~3xlO"''*ergs ' cm"') of ILorenzetti et al 1 (.2009.) spectra, 
taken between 2007 and 2008. 



http://physics.nist.gov/PhysRefData/ASD/lines_form.html 



All lines are spectrally resolved, with FWHM values rang- 
ing from ~100 to ~250kms ', showing mostly broad single- 
peaked profiles. All Brackett line peaks are systematically blue- 
shifted by — 20(+15)kms"', whereas v^ipeak) of the atomic 
lines (blended and low S/N lines excluded) is around 0km s', 
within a +15kms"' uncertainty. On the other hand, the bright 
Pa/S and Pa5 lines are double peaked, with strong red-shifted (at 
~30kms"') and weak blue-shifted (at — 270 km s"') peaks. The 
[Fe n] lines have both high- and low- velocity components (H VC 
and LVC, respectively) in both red- and blue-shifted jet lobes 
(i. e., the bright iron lines show a peculiar four-peaked profile). 
Values of the [Fe n] peak velocities, averaged over all the de- 
tected lines, are -265 and -99kms ' (blue-shifted HVC and 
LVC, respectively), and 165 and 50kms ' (red-shifted HVC 
and LVC, respectively). The [Fen] emissions of the HVCs are 
spatailly extended in our spectral images, indicating that the 
blue- and red-shifted lobes point northwards and southwards, 
respec tively, as also observed in the large-scale CO outflow 
maps (lArce & Goodmanl2002bl) . 
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Fig. 2. LUCIFER z- 
labelled. 



J-, H- and .ff-band MR spectrum of PV Cep normalised to the continuum. The most prominent features are 



4. Analysis and Physical Properties of PVCep 

4.1. Extinction toward tlie source and tliejet 

To obtain an accurate value of the visual extinction (Ay) towards 
PV Cep stellar photosphere and its circumstellar region, we use 
H I and [Fe ii] lines, respectively. The hydrogen atom in the NIR 
traces regions very close to the central source, whereas the iron 
forbidden lines delineate the inner jet from few to several hun- 
dred AUs. For each species, we employ transitions arising from 
the same upper level, namely Bry and Pa(J, and [Fe ii] lines at 
1.257 and 1.644//m. We assume that the emission arises from 
optically thin gas, thus their observed ratios depend only on the 



differential extinction. The theoretical values are derived from 
the Einstein coefficients and frequencies of the transitions. [Fe ii] 
transition probabilities were taken from Nussbaum er & StorevI 
(1988). We adopt the Rieke & L ebofskv (1985) extinction law 
to correct for the differential extinction and compute Ay. From 
the Biy/Pa6 ratio we infer Av=10.8±0.5mag. We then calcu- 
late the [Fell] 1.25/1. 64 jum ratio for both HV and LV com- 
ponents in the red- and blue-shifted lobe of the jet, obtaining 
Av=10.1+0.7mag and Av=10.4±1.9mag (HVC and LVC red), 
Av=6.3±0.5mag and Av=7.1±0.9mag (HVC and LVC blue). 
HV and LV extinction values are identical (within the error bar) 
in each lobe, thus we derive a weighted mean for the extinction 
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of the red-shifted lobe (Av(re(i)=10.1+0.7mag), and the blue- 
shifted lobe (Av(^/Me)=6.5+0.4mag). The difference in the ex- 
tinction between the red- and blue-shifted lobes is due to the cir- 
cumstellar matter (mostly caused by the disc in PV Cep), and it 
is us ually observed in pr otostellar jet s close to the so urce (see 
e. g.. iGarcia Lopez et all 12008; Meln ikovet al.ll2009l) . Finally, 
we note that the Ay towards PV Cep varies with time, due to 
circumstellar extinction variability, related to dust condensation 
in the inner disc region after the outburst and/or to the large 
amount of du st hfted by the je t from the disc during and after 
the outbursts jKun et alj|2011h . Our values are consistent with 
those found in literature, ranging between 9 and 14.5 m ag (e. g., 
iLorenzetti efal. 2009: Connellev et al. 2009: Kun et al.ll201lh . 



4.2. Stellar parameters 

As already mentioned in Sect. 13. II our spectrum does not provide 
any firm conclusion on the spectral type of PV Cep. However, 
there are several clues from our data and from the literature in- 
dicating that the central object is likely an intermediate-mass 
young s tar (2-3 Mc?,), with a r elatively early spectral type (F-A, 
see also ILorenzetti et al.l201 iL and references therein). In partic- 
ular, it is worth noting: a) the presence of H2O MASER emis- 
sion jTor relles et al. 1986 ; Marvel 2005) with a high isotropic 
luminosity (~3x10~^Liti: iTorrelles et aP 11986). which implies 
that t he exciting source is, at least, an intermediate-mass star (see 
e. g.. lAnglada et alj|1996l): b) the presenc e of a massive circum- 
stellar disc (~0.8 MrT\ 'Hami douchel20 1 Ol) . uncommonly massive 
even for a young Herbig Ae star. 

Additionally, there are at least two other relevant indica- 
tions in our spectrum supporting that the central object is an 
intermediate-mass young star: the presence of UV pumped lines 
in emission, and the inferred jet velocity. It is worth noting that 
UV pumped lines have been also observed in low-mass erup- 
tive stars. However, they were d etected only in the outburst 
phase (e.g. EX Lup, SpT M; see iKospal e"tai]|201 lb*). In the 
spectrum of PV Cep, these features are observed during a lower 
accretion state (likely close to the quiescent phase), supporting 
the idea that they are excited by the stellar photosphere. An effec- 
tive temperature (r„) of 8 000 K (with density Nh of lO** cm"^; 
iKun et al]|2011h is needed to produce the observed Fen UV 
pumped lines (see e. g., Joha nsson & Letokhov 2004 , 2007, and 
references therein). Therefore, it is quite likely that Ta is the tem- 
perature of the stellar photosphere (Tefr), which excites the Fe 11 
emission in the circumstellar region close to the disc. 

The high jet velocity (~600kms"', see Sect. 14.4b implies 
that the central source is massive enough (2-3 Mq) to accel- 
erate the ejecta. For example, assuming that the jet velocity 
Vjei is close to the escape velocity, i.e. d,> / ~ (2 GM,/j?,)'^^, 
and M»=2.6M0, R, ^2.9Rq (Th e et all [1994) . we obtain 
Vjet ~580kms"', very close to the value measured in Sect. 14.41 

We therefore favour the hypothesis of PV Cep being an em- 
bedded and young Herbig Ae star (SpT A5; Cohen et al. 1981). 
The adopted stellar parameters along with their references are 
summarised in Table [T] 



4.3. Accretion properties 

Several empirical relationships have been proven effective in de- 
riving YSO accretion luminosity from the dereddened line fluxes 
of different accret ion tracers (see, e. g. Muzerolle et al. 1998a, 
l2001l:ICalvetetal.[ 2004: Natta et al. 2006). We use Pa^, Bry line 
dereddened fluxes (adopting Av=10.8mag, see Sect. 14.11 1 and 



Table 1. PV Cephei adopted stellar and disc parameters 



Stellar Parameter Value Reference 

Distance 

M. 
R, 
SpT 

T, 

Mdisc 



500 pc 


1 


2.6 M© 


2 


2.9 Rq 


2 


A5 


1 


8300K 


1 


0.8 M© 


3 


62°+4° 


3 



Re ferences. (1) ICohen et al.l dl981l) : 
rSl lHamidouciia iiom 



(2) iThe et al.l ( fT994b : 



the fo llowing empirical relationships from ICalvet et aP (l2000l 
I2OO4I) . which have been specifically c ahbrated for protostars of 
intermediate mass JCalvet et al.ll2004l) : 



LogiLacclLo) = 1.03 X Log{LpaplLQ) + 2.80 
Log{L„,JLQ) = 0.90 x LogilBry/LQ) + 2.90 



(1) 

(2) 



As a result, we obtain Lacc(.Pap)-^-^LQ, and 
Lacc{Bry)-5.6LQ. For our analysis, we therefore assume 
an average Lgcc value of 5.7 L©. We note that using the rela- 
tionships from iMuzerolle et al.l (Il998al) . which are calibrated 
on solar-mass YSOs, we would get larger values which are 
not self-consistent (namely, Lacc~^ and 26 L©, from the Pa^S 
and Bry lines, respectively). Moreover, it is worth noting 
that Lorenzetti et al. (2009) derived LgcdPap ) and LgcdBry) 
from t heir multi-epoch spectroscopy using IMuzerolle et al.l 
(Il998ah relationships. The resulting LacAPaP) and Lacc(Bry) 
values systematically diverge by a factor of three to five (see 
their Tab. 12), too large even considering the uncertainties on 
the measured fluxes. On the other hand, their data become 
self-consistent (within the error bars) if Eq. [T]and|2]are adopted. 
In this case their va lues are also coincident with Lace estimates 
of iKun et alJ ( l201lh . ranging from ~80L©, during the outburst 
in 2004 down to ~40L© in 2010. These values are up to one 
order of magnitude larger than ours, indicating that Lgcc of 
PV Cep is still decreasing. 

Finally, we get an estimate of the mass accretion rate by com- 
bining the obtained accretion luminosity and the adopted stel- 
lar parameters. Given that Lace is the energy released by mat- 
ter ac creting from the co-r otational radius (at ~5 R,) onto the 
YSO (iGuUbring et al.l[T998l) : 



.c~GM,Macc{l-RJRin)/R, 

Therefore, M^cc is given by: 



M„, 



1.25RJGM, 



(3) 



(4) 



The resulting value is 2x10"^ M© yr"', i. e. at least one order 
of magnitu de lower th a n wha t was measured bv ILorenzetti et all 
(1200 9) and 'K un et all (1201 ll) during and immediately after the 
2004 outburst. This indicates a decrease in the PVCep mass 
accretion rate since its maximum in 2004-2005. The current 
value is, nevertheless, still about one order of magnitude higher 
than the typical v alues of pre-main sequence (PMS) stars o f 
similar mass (e. g.. lCalvet et al.'2004: ' Garcia Lopez et al.l2006h . 
This implies that accretion at high rates is still ongoing, as also 
pointed out by the presence of the C O lines in emission (eg., 
iDavis et alUlQllllBiscava et al.lll997h . 
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Fig. 3. Position-velocity diagram of PVCep around the 
1.644yL(m [Fell] line in the H-band. The continuum has been 
subtracted. The radial velocity on the X-axis refers to the [Fe ii] 
line and it is computed in the LSR. [Fe ii] blue- and red-shifted 
high- and low-velocity components as well as the Br 12 posi- 
tion are marked. The distance from the source is indicated on 
the Y-axis. The YSO position was derived through a gaussian 
fit to the continuum before its removal from the spectral image. 
Positive distance is computed northwards (along the blue-shifted 
lobe direction). The contour levels of the spectral image are 3, 
10, 20, 30, 40, 50, and 60x the standard deviation to the mean 
background. 



4.4. Jet kinematics and dynamics 

Figure [3] shows a continuum-subtracted spectral-image of 
PVCep around the [Fen] 1.644yL(m line in the H-band. Four 
peaks are detected: two blue-shifted (HVC-blue and LVC-blue) 
and two red-shifted (HVC-red and LVC-red). Red- and blue- 
shifted LV components are observed only in the brightest lines 
(at 1.257 and 1.644/im). They are not spatially resolved, ob- 
served on source within ~0.8", or ~400 AU assuming a distance 
of 500 pc. On the other hand, the HVC emission extends up to 
~2'.'5 away from the source, tracing part of the jet close to the 
YSO (the wiggling geometry of the jet did not allow us to en- 
compass the jet full extension). The HVC red-shifted emission 
peaks at ~l" (or ~500 AU) from the source, whereas we do not 
detect any obvious knot emission along the blue-shifted part of 
the jet, although the emission is clearly elongated. Indeed the 
actual jet PA. slightly differs from the adopted slit PA., and 
our slit encompassed the southern jet lobe better t han the north- 
em lob e. This can be clearly seen in Figure 5 of iHamidouchel 
(l2010h . where blue- and red-shifted outflow lobes unequivo- 
cally show different position angles. [Feii] profiles of the 1.257 
and 1.644//m lines, normalised to the continuum, are shown in 
Figure |4] HVC profiles are well defined, with average peak ra- 
dial velocities of — 265 and 165 km s"' and an average decon- 
volved line widtqj (An) of ~70kms '. On the other hand, the 
LVC profiles are less delineated, with average peak radial veloc- 
ities of — 100 and 50 km s"', and an average Ay of ~ 160 km s"'. 
Figure |4] indicates that both HVC and LVC peak velocities are 
asymmetric (i. e. the blue-shifted peak velocities are larger than 
the red-shifted counterparts) by a factor of ~ 1.6 and ~2, respec- 
tively. Moreover, as already seen in Sect. 14.11 the decreasing in- 
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Fig. 4. [Fell] 1.257 and 1.644yum line profiles, normalised to 
the continuum, plotted on a velocity scale in the LSR. The 
Br 12 line has been fitted with a Gaussian profile and removed 
from the 1.644/im line profile. Average peak velocities at -265 
and -99kms"' (blue-shifted HVC and LVC, respectively), and 
165 and 50kms"' (red-shifted HVC and LVC, respectively) are 
marked. 



tensity ratio between the two line profiles implies that the visual 
extinction increases from the blue- to the red-shifted lobe. 

Interferometric observations of 'H amid ouchel (1201 0*) provide 
us with a good estimate of the PV Cep disc inclination (62°+4°), 
indicating that the system axis is relatively close to the plane of 
the sky (28°+4°). From the radial velocity, we can thus derive the 
tangential and total velocity of the jet (iitg and iitot, respectively) 
as well as estimate when the observed knot was ejected. Inferred 
tangential velocities are 500 and 310kms"' for the HV compo- 
nents (blue and red lobe, respectively), and 190 and 90kms ' 
for the LV components (blue and red lobe, respectively), which 
translate into iitot of 565, 350, 210, and llOkms"'. Thus, if we 
assume a iitg of 3 10 km s' for the HVC of the red-shifted jet, 
we infer that the observed knot has been ejected between 7 and 
8 years before our observations (i.e. between 2004 and 2005). 
Therefore, the observed emission is likely the ou tcome of the 
2004 outburst (lElek & KunllJoTol: iKun et atllloTlh . 

4.5. Jet piiysical parameters 

The various [Fe ii] lines also allow us to derive the physical prop- 
erties of the jet close to the source. These lines come from tran- 
sitions among the a^D, a^F, a^D, and a^P terms (see Column 7 
of Tab. Q)- Transitions from the first three levels have similar 
excitation energies {E'^ ~1 1 000-12000 K), but they have diffe- 
rent critical densities («cr ^lO'^-lO^ cm"-'). Therefore, their line 
rati os can be used to derive the gas electron density (n^) (e.g., 
Nis ini et alj|2002t iTakamTet al. 2006). On the other hand, tran- 
sitions originating from the d^P term have a higher excitation 
energy (~20000K), and they can thus provide an estimate of 
the electron temperature (Te), when combined with lines from 
different energetic levels. 
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To infer T^ and n ^ of the jet, we use a non -LTE 
model ("Nisi ni et alJlIool iGarcia Lopez etalJl2008l 12010^ that 
considers the first 16 fine-structure levels of [Fen] and thus 
also includes the lines presented in this paper. Our model em- 
ploys transition probabilities from lNussbaumer & Storevl(ll988l) 
as well as level energies and rate co efficients for electron colli- 
sions from lZhang & Pradhanl(ll995l) . 

We first construct a diagnostic diagram (top panel of 
Figure |5]l using line ratios sensitive to electron temperature 
(0.927 //m/1. 257 jum, X axis) and density (1.600 //m/1.644yum, 
Y axis). The logarithmic theoretical ratios are plotted as a func- 
tion of Te between 10000 and 20 000 K (red soHd lines) and 
He between 10000 to 60000cm"^ (blue dotted lines). We then 
plot the logarithmic dereddened line ratios observed in the blue- 
and red-shifted HVC of the jet (triangular and squared symbols, 
respectively), using the Ay values found in Sect. 14.11 Because 
we do not detect any red-shifted emission in the 0.927jum line 
(due to the high visual extinction), we assume the same elec- 
tron temperature for both lobes. However, this assumption is not 
necessarily correct, because the two lobes have different veloc- 
ities (Sect. 14.41 ). and, possibly, different excitation conditions. 



The derived T^ and «e values for the blue lobe are 17 SOO^^qqIJ K 
and 45 OOO^^jjJJqq cm"-', whereas a slightly lower density value of 
30 000+^5000 cm-3 is inferred for the red lobe. 

A further estimate of the n^ in both lobes is obtained from 
the 1.677yum and 1.644/vm line ratio (see Figure |5] bottom 
panel). For re=17 500K, we derive n^{blue) between 18 000 
and 32000cm-^ and n^ired) between 11 000 and HOOOcm^^ 
Although lower, these estimates are consistent with those previ- 
ously obtained, confirming that the electron density of the red 
lobe is between 1 .5 and 2 times lower than the value obtained 
in the blue part. Thus we get a weighted mean for the electron 
density of 30000cm"^ and 15 000cm"^ (blue and red lobe, re- 
spectively). 



4.6. Mass ejection rate 

Jet physical and kinematical parameters previously inferred per- 
mit us to evaluate the mass ejection rate (Mjet) in both lobes. 
Moreover, by means of our velocity-resolved observations, we 
can infer Mjet in the different velocity components, probing the 
structure of the jet/wind itself. As the [Fe ii] emission is optically 
thin, its luminosity provides us with an estimate of the total mass 
of the emitting gas (Mgas). The mass ejection rate will be then 
given by MgasVtIk, where Vi is the tangential velocity and k is 
the jet extension on the sky plane. The relation between Mgas 
and [Fell] luminosity is Mgas - t^mHinnVgas), where ju = 1.24 
is the average atomic weight, inn and nu are the proton mass 
and the total density, and V is the volume of the emitting re- 
gion. The term {riHVgas) can be reformulated as the luminos- 
ity of a particular transition (the 1.644/im line, in this case), or 
nnVgas = Li,M4(hvAifiS^^-^y\ where L1.644 is the luminos- 
ity of the 1.644yum line, for the selected transition, Ai and fi are 
the radiative rate and the fractional population of the upper level 
of the transition, y is the ionisation fraction of the considered 
species with a total abundance of ^p with respect to the hy- 
drogen. To obtain the luminosity of the 1.644yum line, we first 
deredden the fluxes of the four observed components, accord- 
ing to the Ay values found in Sect. 14.11 The fractional popula- 
tions were then computed using a constant value of 17 500 K for 
Te and Wg values of 30000cm"^ and 15 000cm"^, for the blue 
and red lobes, respectively (see Sect. 14. 5t . We also assume the 
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Fig. 5. Electron temperature and density estimates. Top panel: 
Diagnostic diagram employing line ratios sensitive to elec- 
tron temperature (0.927/im/1.257jum, X axis) and density 
(1.600jum/1.644yum, Y axis). The grid displays a Tg range from 
10 000 to 20 000 K (red solid lines), and an «e range from 10 000 
to 60000cm"^ (blue dotted lines). Triangle and square indi- 
cate the logarithmic dereddened line ratios and uncertainties ob- 
served in the blue- and red-shifted HVC of the jet, respectively. 
Bottom panel: Predicted 1.677//m/1.644//m [Feii] hne ratio as 
a function of the electron density We- The solid curve corresponds 
to re=17 500K. The dashed horizontal lines show the range of 
ratios observed for the blue and the red lobe. The intersections 
between the curve and the dashed lines illustrate the range of rif. 
allowed by the ratios observed in the two lobes. 



same Wg values for the LV and HV components. Such assump- 
tion implies that both components have similar physical proper- 
ties, which is unlikely because they have diflFerent spatial exten- 
sion and might originate from two different c ircumstellar regions 
(e. g., a disc wind and a jet). For example, iGarcia Lopez et all 
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(120081 1201 Ol) find that n^ is higher in the LVC of jets from 
Class I YSOs (about a factor of two with respect to the HVC val- 
ues), whereas_thej3pposite is found in more e volved CTTSs (see 
e. g. jBacciotti et a l. 2000; C offey et alJ2008l) . In the former case 
(HeiLVC) > ne(HVC)) a lower Mjet(LVC) value (less then a fac- 
tor of two) is obtained, whereas we get the opposite in the latter 
case (n^iLVC) < n^{HVC)). 

Finally, we assume that all iron is ionised (i. e. ^-\), with 

a solar abundance of 2.8x10"^ dAsplund et al.ll2005l i. e., no dust 
depletion). 



As a result, we obtain Mjet(i/yCi/„(,)=1.5xlO M0yr" 



and MaHVCredY- 



locity components 
tively, whereas we 
M,e,(LyC.„,)=4xlO- 



Jjet 

.jet(LyC„<;)=4xlO-^Moyr-i 
nents of the blue and red lobes 



1.2xlO"^Moyr"' from 
of the blue and red 
get Mie 



ve- 



the high 
lobes, re spec - 
(LyCw„,)=5xlO-'*Moyr-', and 
from the low velocity compo- 
Uncertainties can be estimated up 
to a factor of two. Despite the asymmetry in velocity, both lobes 
have, within the uncertainties, the same mass ejection rates in 
both HV and LV components, whereas the larger Mjet is carried 
by the HV component, as previously observ ed in other jets from 
young and embedded stellar objects (e. g.. iGarcia Lopez et alj 
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Fig. 6. Close-up around the CO continuum-subtracted spectrum 
(in black) with the best-fitting model superimposed in red. (for a 
detailed description of the modelling, see Sect. 14.7b 



4. 7. CO bandhead modelling 

We investigate whether the CO emission is likely to originate 
from a small-scale gaseous disc, interior to the disc structure 
already resolved by iHamidouchd (l2010l) . To evaluate such hy- 
pothesis and constrain the physical properties of the CO emit- 
ting region, we simultaneously model the CO y = 2 - and 
y = 3 - 1 bandhead emission. A detailed description of our 
modelling technique is given in Wheelwright et al. (2010) and in 
lllee et alJ(l2013h . In brief, a geometrically thin flat disc is consid- 
ered, where the excitation temperature and CO surface number 
density decrease as power laws with increasing radius (r) as: 






(5) 
(6) 



where Ti and A^; are the excitation temperature and the CO 
surface density at the inner edge of the disc (/?;). 

The optical depth is given by the absorption coefficient per 
CO molecule x the CO column density (i. e., A^(r) because of the 
thin disc). Our model takes into account up to 7=100 rotational 
levels in local thermal equilibrium (LTE) for the two transitions, 
where the CO energy levels a nd the Ei nstein coefficients ar e 
taken from Farrenq et al. (1991) and from Chandra et al1(ll996h . 
respectively. The disc is split into radial and azimuthal cells. The 
spectrum of each is calculated individually and then summed to- 
gether to form the spectrum for the entire disc. 

We model the kinematics of the emitting region by assum- 
ing Keplerian rotation (uco = y/GMt/r). We allow the disc to 
be inclined and calculate the line of sight velocity for each disc 
cell. The final emission feature is given by summation of the 
individual cell spectra, which is then convolved to the resolu- 
tion of the observations. The best-fit model parameters are de- 
termined by comparing the observations and model using the 
downhill simplex algorithm implemented in idl as the amoeba 
routine. We explored several different initial positions to avoid 
local minima. The stellar mass and radius are fixed parameters 
(as reported in Table [TJ, whereas the Ti, N;, i, inner radius (Ri), 
intrinsic line width (Au) and the p and q exponents are free pa- 
rameters. Model fits are then related to the data employing the 



reduced chi-squared statistic. The uncertainty in the data is de- 
rived from the standard deviation of the flux in the pre-bandhead 
portion of the spectrum. 

Figure |6] shows a close-up of the spectrum around the CO 
bandheads (in black) with the best-fitting model superimposed 
in red (xj-5.2), whilst the fitted parameters are presented in 
Tabled 

Table 2. Parameters derived from the best-fit model to the CO 
bandheads (x^.=5.2). The outer disc radius (Rom) is defined at the 
point in the disc in which the temperature drops below lOOOK, 
thus no error is presented. 



Physical Parameter 



Value 



Rout 

Ti 
i 

Av 
P 

q 



6x10^1=^1 cm-2 

0.3+0.1 AU 

2.9 AU 

5r 



2930+90 

3+25 



7;^^kms 
-0 5+°-4 

■2 Q+1.8 



The resulting parameters indicate that the CO bandhead 
emission can be reproduced using a disc model that has the 
same inclination as the outer disc. Thus the observed CO emis- 
sion should not originate from the outflow. This supports the 
existence of a small scale gaseous disc structure, that extends 
from ~0.3+0.1AU up to ~2.9AU from the source and has 
an inner excitation temperature of ~3000K. These quantities 
are larger than those found in other EXors with lower M,, 
e.g., EXLupi (R i -0.1 AU, Tgo -25 00 K, M, -0.6 M©; see, 
lAspin et al.l2O10HK6spal et alJlJO l IbV or in CTTs (Ri <0. 1 AU, 
^coring ~0.2 AU ), and more simila r to values found in Herbig 
Ae stars (see e. g.. lSalvk et al.ll201 lb . underlining once more the 
uncommon EXor nature of PV Cep. Our CO modelling is con- 
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sistent with a flared irradiated disc, but it is slightly too hot (sev- 
eral hundreds Kelvin). This might be an aftermath of the recent 
outburst in 2004, which heated and remodelled the inner disc 
region. Finally, assumi ng a dust sublimation temperature (Tsub) 
of 1500K (see e.g., Salvket alj|2011.) . from Eq. |5]we can es- 
timate that the sublimation radius (^sub) is positioned at 0.3- 
0.5 AU from the source. The same result (Rsub ~0.4AU) can 
be obtained from Equation 2 in Salyk et al. ( 201ljj assuming 
a sou rce luminosity (L, + Lace) of ~100Lo (see. lLorenzetti et all 
l2009t and references therein). The inferred Rsub v alue is very 
simila r to those typical for Herbig Ae stars (see e. g.. lSalvk et alj 
1201 1). and it is comparable to the value inf erred in EX Lup dur- 
ing flie outburst (-0.2-0.3 AU: lK6spal et al. 201 lb) . 



(M^ii„d((D) - (dr~* , where w - IsinO is the distance of the point 1,0 
from the rotation axis), which 'distributes' the mass load among 
the streamlines, is 5. The electron temperature is constant and 
equal to 10000 K, the mass loss rate is 3x10"^ Moyr"', and 
/3-4, where the parameter /J is a power index in the velocity law 
(ii(/) = uo + (l>oo - V()){1 - h/l'f, where dq and Hoo are the initial 
and terminal values of the radial velocity). A scaling factor (/) 
between the radial and tangential velocities decreases with the 
distance from the star, due the magnetic field intensity decrease 
along the disc. Table [3] lists the main parameters of our best- 
fit model. In Figure |7] we show the Pa/3 and Bry profiles (red 
dashed lines) resulting from our modelling superimposed over 
the observed profiles (black solid lines). 



4.8. HI emission modeiiing 

As mentioned in Sect. |3] the observed lines from the Paschen 
and Brackett series show different shapes and peak veloci- 
ties: i. e., Brackett lines are blue-shifted and single-peaked (at 
— 20kms"'), whereas Paschen lines are double peaked (see 
Figure |7]i, with bright red-shifted maxima (at ~30kms"') and 
weak blue-shifted peaks (at — 270 km s"'). We aim at modelling 
both strengths and profiles of the Paschen and Brackett lines to 
constrain the physical conditions and study their origin. 

In highly reddened objects, the observed emission is the sum 
of direct and scattered light, both dimmed by the circumstellar 
dust. Depending on the geometry of the scattering and emitting 
regions, the line profile of the scattered radiation can strongly 
differ from the profile given by the direct radiation (Grinin et al. 
l2012h . This circumstance has to be taken into account when 
analysing emission spectra of objects like PV Cep. Moreover, 
extinction diff'erentially affects Paschen and Brackett lines ob- 
served in the infrared region of the spectrum. Finally, scattering 
on dust also takes place in a diff'erent way for lines of the two 
considered series (e. g., the scattering coefficient for the Pa/3 line 
is a factor of 2.6 greater than the Bry line one). 

To model the H i lines, we therefore consider direct radia- 
tion, originating from both magnetospheric accretion and disc 
wind, as well as emission scattered by the circumstellar disc. 
Additionally, our model also takes into account a screening ef- 
fect from the disc. An opaque disc can be a screen that shifts the 
emission lines towards the blue part of the spectrum. We assume 
that the disc is opaque beyond the sublimation radius (~30Rt 
for PVCep). Therefore, the remote red-shifted part of the 
wind is partially hidden from the observer A detailed descrip- 
tion of the disc wind, scattering, and magnetospheric accretion 
models employe d here can be found in Grinin & Tambovtseva 
(2011) ( also see IWeigelt et al.ll201 ll) . iGrinin et al.l (12012 ). and 
iTambovtseva et al.l (1200 Ih . respectively, whilst a complete de- 
scription of the full model, which groups the aforementioned 
models, is given in Tambovtseva et al. (in preparation). For our 
modelling, we consider an inclination of 60° (angle between the 
viewing direction and the system axis), as derived from the ob- 
servations (see Tab. [1] and Sect. 14.71 1. We let vary the mass ac- 
cretion and ejection rate values one order of magnitude around 
those found in Sections [4. 3l and 14.61 

The main parameters o f our disc wind mo del are the fol- 
lowing (see also Fig. A.l in lWeigelt et al.ll201 iL for a graphical 
representation of the disc wind model along with the listed pa- 
rameters). The disc wind launching region (w) extends from 3 
to 37/;, (0.04-0.5 AU); half-opening angles (0) of the disc wind 
for the first and last streamlines are 30°-80°; the parameter y 



Table 3. Parameters of the best-fit model to the Pa^S and Bry 
lines. 



Physical Parameter 



Value 




J 
P 

i 



3-37/;, (0.04-0.5 AU) 

30°-80° 

5 

4 

10 000 K 

60° 

3xlO-**M0yr-i 



(L,„„,Jl67t(rRl,yi^ 



Notes. In our modelling, T^ and ; are fixed parameters, whereas the 
remaining parameters are variable. 



We then also check for the influence of the magnetosphere 
on the line profiles appl ying the magnetospheric model from 
ITam bovtseva et all (1200 ll) . and using an Mom/Mace ratio ranging 
from 0.1 to 0.3. 

From the line profile fitting we note that the contribution of 
magnetospheric accretion to the line profile is small in compari- 
son to the contribution of the disc wind. In our modelling magne- 
tospheric accretion slightly modifies the width of the line, which 
is increased less then 10% of its FWHM. Therefore, it is not 
possible to derive a precise estimate of the mass accretion rate 
from the line modelling. On the other hand, both screening effect 
and the scattering shift the line profiles towards the blue side. It 
is worth noting that both Brackett and Paschen lines would be 
then blue-shifted, but the shift of the Pa/S line is partly compen- 
sated by the self-absorption effect from the outflow, because the 
optical depth of this line is larger than that of the Bry line. 

Finally, the additional blue-shifted peaks in the Paschen 
lines could not be reproduced by our model. Possibly, this 
weak emission originates f rom the jet, which is not included 
in our model. Indeed, Whel an et al.l (2004) observed extended 
Paschen emission associated to jet tracers in several CTTs. The 
velocity peaks of the secondary maxima in our Paschen lines 
( — 270 km s"') have the same values of the [Fen] blue-shifted 
HVC. Nonetheless, no obvious red-shifted component is visible 
in the Paschen profiles. It might be overwhelmed by the main 
emission, assuming that the red-shifted emission of the jet is also 
asymmetric. This would explain the faint bump observed on the 
red-shifted side of the line (see Figure Q left panel). Otherwise, 
its non-detection might be due to the opaque disc, which blocks 
the emission from the counter jet. In this case, the Paschen emis- 
sion from the jet should be very close to the source. Assuming 
that the disc is opaque up to lOOAU, it should be closer than 
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Fig. 7. Pa/3 and Bry observed (in black) and modelled (red 
dashed lines) profiles (see Sect. 14.8b . Lines are normalised to 
the continuum. 

~60AU (assuming an inclination angle of 30° with respect to 
the plane of the sky). 



5. Discussion 

5.1. The accretion/ejection connection 

The detection of an emerging knot, which can be clearly linked 
to the PVCep outburst event in 2004, is indeed an interest- 
ing result, which supports the hypothesis that enhanced accre- 
tion also generates augmented ejection along YSO outflows. 
Boosted winds after outbursts were alr eady observed in Z 
CMa (Benistyetal. 2010) , V I 647 O ri (Brittain etal] |2010|), 
and EX Lup (Goto et al...2()lTt iKospalet al.. .201ldrbv means 
of H I and CO emission lines. Their intensity has been found to 
be closely related with the accretion rate variability. However, 
it is extremely rare that knot formation is observed along the 
jet immediately after an outburst event from an EXor, be- 
cau se EXors usua l ly do not have jets nor show jet line trac- 
ers. IWhelan et al.l (l2010l) report the first case observed in the 
ZCMa companion, which displayed four 'EXor-like' outbursts 
in the past twenty years. The authors present AO imaging in 
[Fen], showing four knots, whose dynamical ages are compat- 
ible with the timing of the outbursts. Thus our data strengthen 
the idea that knots along the flows may have been produced by 
episodic accretion, possibly similar to the 'EXor-like' outbursts, 
i. e. bur sts of short duration, and not as powerful as those from 
FUors (lloannidis & Froebrichll2012l) . In principle, FUor bursts 
might be so powerful that they disrupt the circumstellar envi- 
ronment and the small-scale magnetic field. Thus the coUimated 
jet/outflow would temporarily switch to a wide-angled wind, 
failing in developing knots along the flow. This might explain 
why no knots have been associated to FUor burst episodes. 

A fundamental quantity in the star formation process is also 
given by the mass ejection to accretion rate ratio (Moui/Macc), 
which provides us with the efficiency of the stellar accre- 
tion process. At variance with the single terms of the ra- 
tio, which, on average, depends on M, and decreases as th e 
source evolves ('Natta etal.ll2006t ICaratti o Garatti et al.ll20T2l) . 
Moui/Macc seems to be relatively constant in time (between 
-0.01 and ~0.1), according to both M HD models and observa- 



tions (e. g.. ICalvell2004t iPudritz et al...2007.) . It has also been 
stressed (e. g., ICabritI 120091 and references therein) that such a 
ratio should not exceed a value of about 0.1, in order to agree 
with the main MHD launching models. Indeed, the measured ra- 
tios are strongly affected by the uncertainty on both Mq,,, and 
Mace estimates, usually up to one order of magnitude. From our 



estimates in Sect. l4.3l and l4.6l we would obtain a ratio very close 
to one, which is too high, even when uncertainties are taken into 
account. However, we should keep in mind that the inferred es- 
timates might not trace two simultaneous events. In particular, 
due to the different sizes of the emitting regions, the different 
cooling times of the two considered species, and the low spa- 
tial resolution of our spectral-images (0'.'8~400 AU), the ejection 
rate derived from the [Fe ii] emission is averaged over ~6 years 
around the outburst event, whereas the H i emission is mostly 
tracing the present accretion rate. Therefore we should instead 
compare Mo„, with an average value of Mgcc during a n d after the 
outburst (few 10 ''Moyr '; 'Lorenz etti et al.l 120091: iKun et all 
|2011). The resulting ratio is then between 0.05 and 0. 1 . A similar 
result is obtained comparing the actual MomlMacc ratio, derived 
from the H i in Sect. |43] and Sect.gjl 



5.2. An ionised jet/wind 

The fading phase of PV Cep gives us a unique opportunity to 
study the jet and outflow relatively close to the source, also pro- 
viding us with indications of the accretion/ejection interplay. 

Our analysis shows strong evidence of both jet and wind in 
the PVCep outflowing material. High-resolution spectroscopy 
of YSO jets usually reveals the existence of two velocity com- 
ponents close to the source, namely a high- and low-velocity 
component. The H VC is associated with the extended collimated 
jet. On the other hand, the LVC is usually confined close to the 
central object (up to 100-200 AU) and it is probably produced 
by a disc wind (see e.g., [G arcia Lopez et al. 2008; Dougadq^ 
120081: ICaratti o Garatti & Eis loffel 2009). Both components ai-e 
usually observed through forbidden emission lines (FELs) in 
CTTSs and in a few H erbig Ae/Be stars (see, e. g., iPvo et al] 
I2OO2I : iDavis et al.ll2003h . whereas molecular hydrogen emission 
lines (MHELs) h ave been also dete cted in less evolved Class I 
objects (see, e. g.. lDavis et alJl200lh . Moreover, the LVC is con- 
fined within ~200 AU in CTTSs, whereas in the younger and less 
evolved Class I objects is usually more extend ed (up to ~1000- 
2000 AU; see iGarcia Lopez et al.ll2008l l2010l) . The former pic- 
ture is more likely in PV Cep, in which the [Fe 11] LV component 
is positioned close to the source (<400AU), whereas the HV 
component is extended. In addition, the LVC (in both lobes) has 
a wide FWHM (Au ~160kms"', see Sect, g^l), which is 2-3 
times larger than the HVC FWHM (Av ~70kms"'). Because 
MHD models predict that the jet terminal velocity is mostly de- 
termined by the Alfven radius, wi thin which the magnetic field 
lines act as lever arms (see e. g., iKonigl & Pudritzl 12000 *). the 
width of the line is clo sely related to the size of the emitting 
region (Pvo et al. 2002). Thus the narrow velocity width of the 
HVC indicates that the line-emitting jet is launched within a nar- 
row region in the disc, whereas the broad velocity width of the 
LVC suggests a wider region of the disc. Moreover, the higher 
peak velocity of the HVC implies that it is launched from a disc 
region more internal than that of the LVC. This indicates that 
the HV and LV components trace the jet and the wind, respec- 
tively. Notably, the described pictu re fits well the CO observa- 
tions in lArce & GoodmanI (l2002ah . Based on their analysis of 
the geometry and kinematics of the molecular emission close to 
the source, the authors predicted the presence of a wide-angle 
wind and a collimated jet. Finally, we also note that the jet is 
fully atomic. Indeed the H2 line at 2.12 fj.m (the only H2 line de- 
tected in the spectrum) is not produced along the jet, which is fast 
enough (~600kms"') to dissociate the H2 and partially ionise 
the medium. At variance with the [Feii] emission, the 1-0 S(3) 
line is not spatially resolved and it is confined within the first 
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400 AU from the source. It shows just one velocity component 
at — 13(±15)km s"', i. e., very close to the rest velocity frame of 
the system, with a narrow FWHM (~7.4 A, i. e. An ~60kms"'), 
pointing to a disc and/or wind origin. The former is more likely, 
due to the lack of a red-shifted component. 

5.3. Origin of the jet/wind asymmetry 

Another interesting and peculiar aspect of the PVCep 
jet/outflow is its asymmetry. Velocities in the southern red- 
shifted and northern blue-shifted lobes are asymmetric in both 
LV and HV components, as well as their inferred electron densi- 
ties. Asymmetric jets from YSOs have already been observed in 
various YSOs (se e e. g., Hirth et al. 1994; Melnikov et al. 2009; 
iPodio et alj 1201 lb . but this is tlie first documented case of an 
asymmetric jet produced during an outburst. 

The origin of such asymmetries still remains an open prob- 
lem. It might be intrinsic to the source and strictly linked to 
the launching mechanism (e. g., asymmetry in the magnetic field 
configuration), or extrinsic, i. e. caused by different physical con- 
ditions in the circumstellar medium, such as an inhomogeneous 
medium (for a detailed discussion see, Matsakos et al. 2012). In 
principle, the former scenario could explain both the observed 
YSO variability (which cannot be triggered by the presence of a 
close companion) and the asymmetries observed in the jet. 

However, the estimated mass loss rate in the two lobes is 
comparable, suggesting that the ejection power is similar on 
both sides of the system, as expected from a magneto-centrifugal 
ejection mechanism. Notably, two previously analysed asym- 
metri c jets (namely, RW Aur and DG Tau; see iMelnikov et al.l 
120091; [Podio et alj|201 1 ) show similar characteristics: velocity 
and density asymmetries and the same Mout in both lobes. The 
authors of both works suggest that the observed asymmetries 
are due to different mass load and/or propagation properties in 
an inhomogeneous environment. The medium around PV Cep is 
highly inhomoge neous, as observe d in the large-scale CO out- 
flow studied by .Arce & GoodmanI (l2002bi) . and partially sug- 
gested by the different Ay values observed along the jet in 
the northern (blue-shifted) and southern (red-shifted) lobes (see 
Sect. 14.1b . PV Cep is located in the northern edge of its dark 
cloud, thus the northern and southern lobes interact with en- 
vironments at very different densities. Hence, the blue-shifted 
ejected matter (northwards) interacts with a much less dense 
medium, producing a highly asymmetric jet and, in turn, a highly 
asymmetric outflow. It is worth noting that the velocity asym- 
metry in the [Feii] LV component implies that such inhomo- 
geneities in the medium density must extend very close to the 
YSO (at least down to 400 AU), and, more likely, down to the 
YSO disc. In this case, the disc would directly interact with 
an inhomogeneous asymmetric ambient medium, producing an 
asymmetry in magnetic lever arms and/or launch radii betwee n 
either sides of the disc jFerreira et al.ll2006l iPodio et al.l 1201 lb . 
Thus, the power generated by the MHD dynamo and conveyed 
by the rotating disc to both sides is identical, but the physical 
conditions of the ejected material may be different on the two 
sides. 



5.4. Circumstellar environment 

Our analysis provides us with an overview of the circumstel- 
lar environment of PV Cep during the pre-quiescent phase, i. e. 
when the outbursting phase is fading and the inner disc region 
is rearranging in a more stable state. The inner region is domi- 




Fig. 8. Sketch (not in scale) of PV Cep's innermost circumstellar 
region. 



nated by hydrogen gas (up to ~0. 14 AU) as delineated by the Bry 
line (see FiglS]). A more extended Hi region (up to few AUs) is 
also traced by the Paschen emission lines. In this region there 
is also evidence of atomic fluorescent emission lines (Fe ii at 1 
and 1.688yum; and Ci at 1.17-1.19//m) and other metallic lines 
(e.g., Fei, Fell, Mgi, Ci, Nai), whose excitation potential (4- 
8eV) is lower than Hi (~13.6eV). According to iKospal et al.l 
( 1201 Ibb . the emitting region is placed where the hydrogen is 
mainly neutral, likely shielded from the shock radiation by the 
circumstellar disc. On the other hand, the fluorescent emission 
indicates that the emitting region is located close to the disc 
surface. Beyond this region, we model a small-scale gaseous 
disc ring, traced by the CO emission, and extending from ~0.2- 
0.4 AU to ~3 AU from the source. The inner region has a temper- 
ature of ~3000 K. This part of the disc is hotter in comparison 
with the typical values found in Herbig Ae stars. This might be a 
result of the heat produced by the outburst. Dust is likely to co- 
exist within this gaseous region, Z?,,,/, being at 0.3-0.5 AU from 
the source. 

Finally, [Fe ii] analysis proves the existence of both wind 
and jet emission close to the source. Our modelling suggests 
that the bulk of the H i emission originates from a wind located 
close to the source. The presence of a secondary maximum in 
the Paschen lines also indicates that a small part of the Paschen 
emission might originate from the blue-shifted part of the jet. 

5.5. How close is PVCep to EXors? 

An interestin g question to answer is how similar PV Cep is to 
other EXors. iLorenzetti et al.l dMll) already pointed out that 
PVCep is not a genuin e EXor according to the classical def- 
inition given by iHerbia (Il989l) . As considered by the authors, 
the main resemblance with the EXor group resides in both spec- 
trum and recurrent outbursts driven by disc accretion episodes. 
On the other hand, the most important discrepancies appear to 
be its higher mass and earlier evolutionary stage with respect 
to other classical EX ors. The remaining differences listed in 
ILorenzetti et al.l(l201 lb (i. e. circumstellar nebulosity, jet/outflow 
presence, massive disc, far-infrared excess in the SED, presence 
of maser and radio continuum emission) can all be a conse- 
quence of the previous two. 

We can now compare our results with the circumsteUar en- 
vironment of the EXor prototype EX Lup (iKospal et al.ll201 Ibt 
iSicilia-Aguilar et al. 2012) to get a more quantitative evaluation. 
Our analysis shows many similarities between the two objects, 
both in the physical structure and the observed spectra. However, 
as already mentioned in Sect. 14.71 and 14.8! the size of the vari- 
ous regions appears to be larger and the excitation conditions 
are more extreme. This is compatible with a hotter, more mas- 
sive object. Another significant difference in PV Cep spectrum is 
the presence of jet tracers. Conical nebulae (namely the outflow 
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cavities), jets, knots and outflows are not present in EXors. This 
indicates that the PV Cep is young and that its circumstellar en- 
vironment has not been cleared by the YSO activity, favouring 
the formation of a jet, shocked knots, and a swe pt out outflow. 

In conclusion, as already pointed out by iLorenzetti et alj 
(1201 ih . PVCep is not a classical EXor because it is not a 
low-mass PMS star It is more massive, young and embedded 
than classical EXors. Nevertheless, it shows the signature of an 
'EXor-like' outburst, which suggests a common origin. We can 
therefore argue that PV Cep is a younger, high-mass counterpart 
of classical EXors. Indeed, the discovery and study of more em- 
bedded outbursting YSOs in the near future will increase our 
knowledge of and statistics on these phenomena and possibly 
confirm that enhanced accretion processes are common to all 
YSOs, despite the difference in their mass or evolutionary stage. 
Therefore, on the basis of the new discoveries, it is quite likely 
that the two main sub-classes (FUors and EXors) will be revised. 

6. Conclusions 

We have presented near-IR medium-resolution spectroscopy of 
PVCep taken with LUCIFER at the LBT in 2012, i.e. eight 
years after its last outburst. By means of the many emission lines 
detected in the spectrum, we are able to trace both accretion and 
ejection activity, deriving the physical and dynamical conditions 
in its circumstellar region. The main results of this work are the 
following: 

- The NIR spectrum of PV Cep displays several strong emis- 
sion lines and a steeply rising continuum without photo- 
spheric absorption features. The most prominent features 
in the PVCep spectrum are circumstellar features, originat- 
ing from YSO accretion or inner winds, as, e. g., H i. Can, 
shocks along the jet ([Fe ii]), inner disc region (e. g.,Na i, CO) 
or chromospheric activity (e.g., Fe i, Fe ii, Mg i, C i). A few 
fluorescent emission lines fromFeii(i. e., at 1 and 1.688/vm) 
and Ci (i. e., 1.17-1.19//m) are pumped by UV photons. 

- Compared to the outburst in 2004 and subsequent photome- 
try, our data show that PV Cep's brightness is fading, indicat- 
ing that the outburst phase is declining. However, the many 
detected lines show that accretion and ejection activities still 
proceed at high rate. 

- The actual mass accretion rate is 2x10"^ M0yr"', i.e. at 
least more than one order of magnitude lower than values 
measured during the 2004 outburst. 

- Among the several emission lines, only the [Fe ii] intensity 
increased since the outburst in 2004. The observed emis- 
sion traces blue- and red-shifted lobes of an asymmetric 
jet/outflow. In each lobe, two velocity components are de- 
tected: a LVC, not spatially resolved (within 400 AU from 
the source), likely tracing a disc wind, and an extended HVC 
tracing the jet. Velocities between the two lobes are asym- 
metric (~1.6 ratio between blue- and red-shifted part). Total 
velocities for the HVC and LVC are -570 and 350km s', 
and -210 and 130kms ', respectively. 

- The observed emission has a dynamical age of 7-8 years, 
indicating that it was produced during the 2004 outburst. 

- Jet physical properties are derived from different [Fe ii] tran- 
sition ratios. The visual extinction decreases moving from 
the red-shifted (Av(reii)=10.1+0.7mag) to the blue-shifted 
lobe (Av(/7/Me)=6.5±0.4mag). We derive an average elec- 
tron temperature of re=17 500K, and electron densities of 
30000cm"^ and 15 000cm"^ for the blue and the red lobe, 
respectively. 



- The measured mass ejection rate is ~L5xlO"^M0yr"' in 
both lobes. These estimates fairly match the high accre- 
tion rate observed during and immediately after the outburst 
(M„„,/M„,, -0.05-0.1). 

- The similar mass loss rate in the two lobes suggests that 
the ejection power is comparable on both sides of the disc, 
as expected from MHD ejection mechanisms. The observed 
asymmetries are thus consistent with an inhomogeneous 
medium, as also indicated by the large-scale CO analysis in 
lArce & Goodmanl(l2002ah . 

- Our modelling of the CO emission hints at a small-scale 
gaseous disc ring, extending from ~0.2-0.4AU to ~3AU 
from the source. This region has an inner temperature of 
~3000K (exponentially decreasing). 

- Our H I line modelling indicates that most of the observed 
emission comes from an expanding wind at re=10000K. 
Brackett and Paschen line profiles are strongly affected by 
scattering, disc screening, and outflow self-absorption. 

- PVCep is not an EXor ob j ect, a ccording to the classical 
definition given by iHerbid (Il989l) . It is more massive and 
younger than typical EXors. Nevertheless, it shows the sig- 
nature of an 'EXor-like' outburst, suggesting a similar out- 
burst mechanism. 

In conclusion, EXors and related objects offer a unique op- 
portunity to investigate the fundamental accretion/ejection 
mechanism, because the physical properties of the accretion 
process and, only few years later, the corresponding ejection 
event can be related and studied in detail. 
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Table .1. Observed emission lines in the spectrum of PV Cephei 



^obs 


F 


AF 


w. 


FWHM 


S/N 




ID (term) 


^vac 


V, {peak) 


(A) 


(10-15 


ergcm-^s-i) 


(A) 


(A) 


ratio 








(A) 


(kms-') 


9173.9 


11.4 


2.1 


-10.4 


6.4 


5 




Fen 




9175.2 


-18 


9230.4 


6.0 


1.0 


-6.2 


5.3 


6 




Pa 9 




9231.6 


-14 


9261.3 


2.3 


0.5 


-2.1 


3.1 


5 


[FeiiKa^Pi/j- 


a'Fy2) 


9270.1 


-262 


9994.4 


1.7 


0.5 


-1.1 


4.3 


3 




Fei? 




9994.0 


38 


10000.3 


3.2 


0.7 


-2.1 


3.8 


5 




Fen 




10000.3 


20 


10042.7 


5.8 


0.9 


-1.6 


7.2 


6 




Pad 




10052.2 


-259 


10052.4 


15.6 


1.0 


-8.7 


7.7 


16 




Pad 




10052.2 


31 


11640.5 


3.2 


0.3 


-1.1 


4.3 


11 




Fei 




11641.5 





11661.6 


2.0 


0.3 


-0.8 


4.2 


7 




Ci 




11662.9 


-8 


11673.6 


2.3 


0.5 


-0.8 


6.6 


5 




Ci + Fei 


7 


11672.8 


44 


11692.2 


4.9 


0.4 


-1.8 


5.0 


14 




Fei 




11693.2 





11750.4 


4.0 


0.4 


-1.4 


6.2 


9 




Ci 




11751.5 


-2 


11756.2 


13.0 


0.5 


-4.4 


7.1 


26 




Ci 




11756.6 


14 


11780.0 


1.7 


0.3 


-0.6 


4.0 


6 




Ci 




11780.8 


5 


11785.4 


5.5 


0.4 


-1.9 


6.1 


13 




Fei 




11786.5 


-3 


11830.2 


13.1 


0.5 


-4.3 


6.5 


29 




Mgi 




11831.4 


-6 


11840.9 


10.7 


0.5 


-3.5 


6.9 


22 




Can 




11842.2 


-9 


11885.1 


14.8 


0.5 


-4.7 


5.7 


32 




Pii + 




11886.1 


-1 


11885.1 














Fei 




11886.1 


-1 


11894.9 


1.5 


0.4 


-0.7 


6.3 


4 




Ci 




11896.2 


-7 


11897.5 


2.1 


0.4 


-1.1 


6.7 


5 




Ci 




11899.0 


-13 


11952.2 


6.9 


0.5 


-2.1 


5.9 


13 




Can 




11953.0 


4 


11975.1 


13.9 


0.5 


-4.4 


5.8 


27 




Fei 




11976.3 


-6 


11986.0 


8.8 


0.5 


-2.8 


5.5 


18 




Caii + 




11987.5 


-13 


11986.0 














Sii 




11987.5 


-13 


11993.5 


6.2 


0.5 


-2.0 


5.3 


11 




Caii + 




11994.8 


-9 


11993.5 














Si I 




11994.8 


-9 


12033.6 


12.1 


0.6 


-4.1 


6.3 


19 




Si I 




12034.8 


-5 


12085.9 


5.2 


0.6 


-1.6 


5.8 


9 




Mgi 




12086.6 


8 


12105.7 


4.2 


0.5 


-1.3 


5.2 


8 




Si I 




12106.9 


-4 


12272.9 


6.4 


0.4 


-1.8 


6.7 


16 




Si I 




12274.1 


-3 


12558.2 


16.6 


0.2 


-4.4 


4.0 


67 


[Feii](a4D7/2- 


a^'D^li) 


12570.2 


-263 


12565.2 


9.9 


0.3 


-3.0 


6.8 


23 


[Feii](fl4D7/2- 


fl''£>9/2) 


12570.2 


-96 


12570.9 


3.7 


0.4 


-1.5 


7.0 


13 


[FeiiKfl^Dv/a- 


fl''D9/2) 


12570.2 


40 


12575.9 


7.0 


0.2 


-1.8 


3.4 


34 


[¥ei\\ {a'^D-j 12 - 


a^D<)i2) 


12570.2 


160 


12616.9 


2.2 


0.4 


-0.6 


6.9 


5 




Ci 




12617.6 


9 


12694.6 


2.0 


0.3 


-0.5 


3.3 


7 




Sii? 
Ah? 




12694.2 
12695.0 


32 
15 


12779.0 


1.9 


0.3 


-0.5 


3.1 


7 


[Feii](fl'*D3/2- 


a^D3/2) 


12791.3 


-263 


12785.9 


1.1 


0.2 


-0.3 


4.3 


5 


[Feii](fl4D3/2- 


a^D3/2) 


12791.3 


-101 


12796.8 


1.7 


0.5 


-0.4 


5.3 


4 


[Feii](fl4D3/2- 


a^Di^2) 


12791.3 


154 


12808.7 


4.4 


0.2 


-1.0 


4.1 


20 




Pa;8 




12821.7 


-279 


12821.7 


113.0 


0.8 


-26.8 


9.4 


134 




Pa;6 




12821.7 


25 


12933.8 


3.4 


0.5 


-1.0 


4.3 


6 


[Feii](fl4D5/2- 


a^D,i2) 


12946.2 


-263 


13125.6 


5.5 


0.6 


-1.7 


5.0 


9 




Ah 




13127.0 


-7 


15441.6 


16.5 


1.1 


-4.0 


12.3 


15 




Br 17 + ■ 


7 


15443.2 


-6 


15559.0 


24.0 


1.0 


-4.5 


11.1 


24 




Br 16 




15560.8 


-9 


15624.3 


3.7 


1.0 


-1 


7.5 


4 




7 








15703.1 


24.0 


1.0 


-4.5 


10.7 


25 




Br 15 




15705.0 


-12 


15742.2 


1.9 


0.5 


-0.6 


5.4 


4 




Mgi 




15745.0 


-29 


15750.2 


10.5 


1.0 


-3.1 


11.5 


10 




Mgi + ■; 


\ 


15753.2 


-32 


15769.5 


18.0 


0.8 


-3.0 


8.6 


23 




Mgi + ■; 


» 


15770.0 


15 


15883.0 


38.0 


1.3 


-7.0 


11.0 


29 




Br 14 




15885.0 


-12 


15889.5 


42.0 


1.3 


-7.0 


11.0 


32 




Mgi 




15890.6 


4 


15961.4 


3.3 


0.7 


-0.6 


5.6 


5 




Sii + ? 




15964.0 


-25 


15965.8 


3.6 


0.7 


-0.6 


5.8 


5 




Fei 




15966.9 


4 


15983.3 


6.3 


0.8 


-1.0 


7.0 


8 


[Feii](fl4D3/2- 


■ a%i2) 


15999.1 


-272 


16007.1 


5.3 


0.7 


-0.9 


5.9 


7 


[Feii](fl4D3/2- 


■ a%i2) 


15999.1 
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-^obs 


F 


AF 


w. 


FWHM 


S/N 




ID (term) 




^uac 


v.ipeak) 


(A) 


(10-15, 


ergcm-^s-i) 


(A) 


(A) 


ratio 








(A) 


(kms-i) 


16111.6 


41.3 


1.4 


-6.7 


11.6 


30 




Br 13 




16113.8 


-16 


16382.3 


5.4 


0.8 


-0.8 


6.5 


7 




Sii + Fei 




16382.7 


15 


16409.3 


55.1 


1.5 


-8.4 


9.3 


30 




Br 12 




16411.7 


-20 


16424.1 


21.1 


0.8 


-3.1 


6.0 


27 


[Fe II] {a%i2 - 


fl'*F9/2) 


16440.0 


-265 


16433.2 


14.6 


0.9 


-2.1 


10.0 


8 


[Fe II] (fl't)7/2 - 


fl'*^9/2) 


16440.0 


-99 


16442.0 


7.5 


0.6 


-1.1 


10.0 


5 


[Fe II] (fl't)7/2 - 


fl'*F9/2) 


16440.0 


61 


16447.7 


13.4 


0.5 


-2.4 


4.5 


28 


[Feii](fl4Dv/2- 


a^^9/2) 


16440.0 


165 


16626.0 


2.4 


0.6 


-0.3 


4.8 


4 


[Feii](a4Di/2- 


a^Fsii) 


16642.2 


-267 


16757.3 


4.8 


0.6 


-0.7 


4.8 


8 


[Feii](a4D5/2- 


a%ii) 


16773.3 


-262 


16781.3 


4.3 


0.8 


-0.7 


7.0 


5 


[Feii] (fl4£)5/2- 


a%/2) 


16773.3 


167 


16808.9 


62.0 


1.4 


-8.4 


11.4 


45 




Br 11 




16811.2 


-16 


16875.8 


6.2 


1.3 


-0.9 


8.5 


5 




Fell 




16877.8 


-11 


16893.4 


18.3 


1.3 


-2.5 


9.0 


15 




Ci 




16895.0 


-4 


17111.7 


25.4 


1.4 


-3.5 


9.8 


19 




Mgi 




17113.3 


-4 


17327.6 


5.4 


1.2 


-0.9 


8.6 


5 




? 








17337.5 


4.4 


1.6 


-0.7 


11.5 


3 




Fell 








17344.5 


6.0 


1.9 


-1.2 


13.5 


3 




Nai 




17345.8 


2 


17364.3 


85.4 


1.8 


-10.6 


12.8 


48 




Brio 




17366.9 


-21 


17406.5 


12.0 


1.2 


-1.7 


8.7 


10 




Fei + ? 




17406.9 


17 


21215.6 


12.6 


1.8 


-1.5 


7.4 


9 




H2(1-0S(1)) 


21218.3 


-13 


21658.1 


193.0 


3.0 


-16.6 


15.5 


66 




Bry 




21661.3 


-20 


22061.0 


25.6 


2.2 


-2.1 


11.2 


11 




Nai 




22062.4 


6 


22087.4 


21.8 


1.6 


-1.8 


8.2 


13 




Nai 




22089.7 


-7 


22711.8 


18.8 


3.6 


-1.5 


16.2 


5 




? 








22937.2 


267.0 


3.6 






5 




CO (2-0) 




22935.0 


53 


23235.0 


200.0 


3.6 


-19.0 


41.0 


5 




CO (3-1) 




23227.0 


128 



